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Abstract 

A closed-form exact analytical solution for the g-state Potts model on a 
ladder 2 x co with arbitrary two-, three-, and four-site interactions in a unit 
cell is presented. Using the obtained solution it is shown that the finite- 
size internal energy equation [6] yields an accurate value of the critical 
temperature for the triangular Potts lattice with three-site interactions in 
alternate triangular faces. It is argued that the above equation is exact 
at least for self-dual models on isotropic lattices. 
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1 Introduction 



The methods aUowing to extract information about a multi-dimensional system 
from solutions of their lower dimensional counterparts play an important role 
in statistical physics. One of the most well-known examples of this kind is the 
finite-size scaling approach [|[ ||. 

There are cases which evoke particular interest when the critical properties 
of a system experiencing a phase transition can be determined exactly from data 
only for its subsystems. For instance, for the Ising strips the intersection point 
of the locus of partition function zeros in a complex temperature plane with the 
real positive axis yields the exact value of the critical temperature for the two- 
dimensional Ising model Exact critical temperatures for the 5 = 1/2 Ising 
models on square, triangular, honeycomb, and centered square (Union Jack) 
anisotropic lattices are obtained by using strip clusters when an effective field 
is applied to one side of a strip only [|j . Another exotic way of estimating the 
critical point of the square- lattice Ising model was proposed in |^ . The authors 
of this paper showed that, in the quasi-diagonal form of a transfer matrix of 
finite width strip, all coefficients of the characteristic equation for the sub-block 
containing the largest eigenvalue have extremum located precisely at the exact 
value of phase transition temperature of infinite lattice. 

In the present paper we concentrate our attention on the method to calculate 
the critical temperature, proposed by Wosiek |^ (see also |0, ||, ||, |lO[ |ll|, p^). 
The author introduced a maximum criterion for the ratio of moments of 
the transfer matrix, and obtained the following equation for determining the 
position of critical point in a d-dimensional system: 

ui{K,) = U2{K,). (1) 

Here ui and U2 are the internal energies of {d~ l)-dimensional and two coupled 
[d — l)-dimensional subsystems, respectively, and Kc is the critical coupling 
(normalized inverse critical temperature) of the d-dimensional system. 

It is remarkable that at d = 2 Eq. (|l|) , shown in Q , yields the exact value of 
Kc for the isotropic square and triangular Ising lattices, as well as for the three- 
site Potts model on the square lattice also with isotropic interactions. Later on 
some other models were added to the list. Now it includes another isotropic 
Baxter model (two square Ising lattices coupled by four-particle interactions), 
the Baxter- Wu model (triangular lattice with three-site interactions of Ising 
spins) [0, and q-state Potts model on the isotropic square lattice with arbitrary 
value of q . Physical nature of Eq. (|l]) may be elucidated when it yields either 
(a) an exact solution, or (b) approximate estimation, or (c) does not give any 
solution at all for a given model. The more such examples are found, the better 
we understand its sense. 

For a two-dimensional system, Eq. (0) connects the internal energies of in- 
finitely long linear and double chains. Therefore, to obtain a possibility to test 
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Eq. rigorously, it is necessary to have analytical solutions for such subsys- 
tems. 

In the second section of the present paper, we give an exact analytical solu- 
tion for the two-chain Potts strip with a large number of independent parame- 
ters. As a special case, it contains a solution for the linear Potts chain. 

Our solution for the double Potts chain enables us to cover all earlier re- 
vealed cases in which Eq. (|^) reproduces exactly the critical temperatures for 
the two-dimensional models of Ising, Baxter- Wu, and Potts. Besides, we dis- 
cover (section ||) a new model for which Eq. (|l|) yields the exact result. This is 
the g-state Potts model on the triangular lattice with pure three-site interactions 
in half of the triangular faces jl^ . 

In the fourth section of this paper a discussion of results is presented. In 
particular, it is shown that duality is a sujjicient condition for the validity of 
Eq. (|l]) for the isotropic spin lattices. 

Finally, in section 5, we summarize the results obtained in the work. 



2 Solution of double g-state Potts chain with 
symmetry 

Consider a two-chain (ladder) lattice with spin variables al attached to its sites 
(i = 1, 2 is the chain index, and 1 = 1,2,3,... labels the sites in the longitudinal 
direction of the ladder); the spin variables take the values 1,2, ... ,q. 
Let us write the Hamiltonian of a system in the form 

n = -J2H{al,af;al„al,). (2) 

The locality of interactions in the above Hamiltonian allows one to introduce 
the transfer matrix V with elements 

(CTi,CT2|y|CTi,(T2) ^ exp[H{ai,a2;(Ji,(T2)/kBT] (3) 

(T is the temperature and fc^ is Boltzmann's constant) and reduce the problem 
of calculating the free energy density / of infinitely long strip to a search for 
the largest eigenvalue Ai of the matrix V: 

/=^lnAi. (4) 

The transfer matrix (^) has the size x q^. It is real and, moreover, all its 
elements are positive, but the matrix is, generally speaking, non-symmetrical 

To solve the eigenvalue problem for the transfer matrix, we shall use a group- 
theoretical approach (in this connection see, e.g., Ref. where such an ap- 
proach was applied to a quasidiagonalization of a transfer matrix of Ising model 
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on parallelepipeds L x L x oo). In order to obtain a solution for the two- leg 
spin ladder (in which we are particularly interested) in a most general form, 
we will go in reverse order. Namely, we first select a symmetry group in the 
space |(Ti,(72), that would enable us to quasidiagonalize the transfer matrix up 
to sub-blocks secular equations which can be solved analytically; only then we 
expand the Hamiltonian density H into a series on invariants of the symmetry 
group. 

Let us take a model which is invariant, for example, under transformations 
of the symmetric group Sq of degree q. For the Potts model this means, as a 
matter of fact, that we are dealing with a system in the zero external field. But 
fortunately the field is not required to test Eq. (||). 

It is known (see, for example, [^) that the largest eigenvalue of a transfer 
matrix is located in the sub-block of the identity irreducible representation (IR). 
In accordance with group theory, the basis vectors ipi of the identity IR can be 
obtained acting by the permutation operators of a group Sq successively on the 
orths |1, 1), |1, 2), . . . , \q, q). Acting by elements of symmetric group first on the 
orth |1, 1) and then on the orth |1, 2), we find that the two linear combinations 
obtained involve all orths. The normalized basis vectors are given by 



(the prime at the second sum indicates that the terms with i — j are omitted) . 
Hence, the sub-block of identity IR has sizes 2 by 2 and, therefore, its eigenvalues 
(one of which is Ai) can easily be obtained by solving an algebraic equation only 
of second degree. Notice that if one takes the group x Cs (Cs is the group of 
mirror refiections in the plane going between the chains of the two-leg ladder), 
the sub-block corresponding to the identity IR will again have the size 2x2 
and therefore this symmetry, which only reduces the number of independent 
parameters in a Hamiltonian, does not justify itself in the given case. 

We now represent the Hamiltonian (^) in the form of a sum of terms which 
are invariant under transformations of the group S^: 






+ Ji5^1„2 1 2 ]. 

"l"l 



(6) 



The Kronecker symbols entering here are defined as follows: 
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The structure of two-site couplings in the Hamiltonian (^) is shown in Fig. 1. 
Matrix elements of the original transfer matrix are written as 

{ai,(J2\V\cj[,(j'^) = exp[Xi(5^^^^ +K25^^^'^ + ]^KQ{5a,a2 +^a[aO 

(TiCT,-, <y2^i <-'l<^2 l7i<T^ (71(7-^ ^'^'-^2 

+k6„,aJ^ • +K 5^ , +K4S„ '5„,aJ^ (8) 

where Ko = Jo/ksT, Ki = Ji/keT, K2 = J2/kBT, K = j' /keT, K" = 
j" /ksT, K3 - Ja/keT, = Jg/fesT, K3 = Ja/keT, = J^keT, and 
K4 = J^/ksT. 

Using expressions (||) and we calculate the matrix elements Qij = 
ipfVipj of the sub-block corresponding to the identity IR: 



hi = [g - 1 -t- exp(is:i +K2 + K + K +Kz + K^ + kz + k^ + K^)] 
>12 = {q- - 2 + exp(ifi + k" + K3) + cxp(i^2 + K + A'3)]e 



1 

Q21 = (g - - 2 + exp(ifi +k' + k'^) + exp(i^2 + K" + K^)]e'^^\ 

Q22 - (g - 2)(g - 3 + + e^^ + e^' + e^") + exp(ifi + K2) + exp(if' + x"). 

(9) 

As a result, we find that the largest eigenvalue of the transfer matrix of the 
double g-state Potts chain with the Hamiltonian (^) reads 

Ai'^ = \{Qii + Q22) + [\{Qii - Q22? + {q- l)Ae'''Y'\ (10) 

where 

A^[q-2 + exp{Ki + K" + ifg) + cxp(X2 + K + A'3)][g - 2 

-Kexp(ii:i +k' + k'^) + eyi^{K2 + k" + K^)]. (11) 

Those variants of the double Potts chains solved earlier |l2], |l^, |l^, Q 
correspond to a particular choice of the interaction constants. Setting Jq = 
J (= J) with all other interaction constants vanishing, we arrive at the solution 
for the linear Potts chain [l9| : 

\^^\K) = e'' + q-l. (12) 



3 Triangular Potts lattice with three-site inter- 
actions on alternate triangle faces 

A large number of independent parameters in the model, solved in the previous 
section, enables us to test Eq. (|l|) for a wide class of two-dimensional spin 
systems. 
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In addition to the cases listed in the Introduction, in which Eq. (|l|) is satisfied 
exactly, let us consider the Potts model on a triangular lattice with three-site 
interactions in each up-triangle (Fig. 2). The position of the critical point in 
this model was found with both three- and two-site interactions p3| . However, 
it is known |Q that for the triangle lattice with pair couplings, Eq. (||) yields 
the exact result only for the Ising case {q = 2). Therefore we discuss the model 
with purely three-site interactions. For this case, 

K, = \n{l + q). (13) 

Let us show that this value satisfies Eq. (|l|) by subsystems in the shape of strips 
with a periodic boundary condition in the transverse direction. 
The internal energy of one-dimensional subsystem is 

u,{K)^^ = [{q-l)e-'' + l]-\ (14) 

Substituting K^^ — with all other interaction constants vanishing, from (^, 
WW, and (O) we obtain an expression for the free energy density of double 



Potts chain: 



f2{K) = \ h\{^^ + q'-l) + [\{e'^ -{q- iff + q{q - l)(2e^ + 
q-2)Y'\ (15) 



The internal energy U2{K) = df2/dK. Differentiating Eq. (g_5|) with respect to 
K one finds an expression for U2{K). 

An analysis shows that the dependences ui{K) and U2{K) have a crossing 
point which lies exactly aX K — K^. = ln(l + q) both for integer and non-integer 
q. The internal energy of a system at critical point is equal to 

uUKc) = ui(ifc) = U2{K,) = i(l + q-1). (16) 

Thus, using solutions only for the linear and double Potts chains, Eq. (|^) has 
enabled us to extract the exact value of Kc for the bulk two-dimensional Potts 
model on a triangular lattice alternating faces which interact by three-site forces. 

4 Discussion 

In the paper ^ Eq. (|l|) was extended up to the form 

ul{K,) = ul'{K,) (i,L' = 1,2,3, ...). (17) 

Here is the internal energy per site of i coupled (d— l)-dimensional subsys- 
tems. In the two-dimensional case L denotes the width of a strip. 
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The validity of condition (^7|) for arbitrary L and L means the absence of a 
"singular" (i.e. depending on L) part of the internal energy density at critical 
point: 

ul{Kc) = const upon L. (18) 

In other words, the amplitudes of all finite-size corrections to the critical internal 
energy of a system Uoo(-f^c) are equal to zero. 

For the square isotropic Ising lattice, the derivative of the inverse correla- 
tion length kl{K) with respect to a temperature-like variable K has a similar 
property MA 



dK 



Ok 



l' 



dK 

K=K^ 



(19) 



i.e. dKL/dK\c does not depend on L. This property has enabled us to determine 
exactly the value of the thermal critical exponent yt {— 1) for this model using 
only the finite-size dataj 20[ pTf 



Equations (|l|) and (|T7) are valid for the ferromagnetic isotropic square Potts 
lattices. These models are self-dual and the critical coupling for them (in the 
anisotropic case) is determined from the condition 

{e^- ~ l){e^y - I) = q. (20) 

For the antiferromagnetic square-lattice Potts model the criticality condition is 

[H 

(e^- + l)(e^« + l) = 4-g, (21) 

where both < and Ky < 0. We performed a verification and found that 
in the antiferromagnetic case the curves ui{K) and U2{K) do not have any 
self-crossing point and therefore Eq. (|l|) does not lead to the exact value which 
follows from Eq. ( pi] ) or even to any approximate estimate for the critical point. 

It is not difficult to show that if a model is self-dual and by this the dual 
point coincides with the original one, then Eqs. @) and (|l^ are valid. 

Indeed, consider for instance the Ising model on the isotropic square lattice 
L X N with toroidal boundary conditions. The partition function of such a 
system has a fundamental property: it is invariant (up to multiplicative factor 
exponentially depending on L) under the duality transformation (see psf): 

ZL.NiK*) = {sin\i2K)-^^'ZL,M{K), (22) 

where K and K* are related by the condition 

tanh/^* = e"2^. (23) 

(We here used another normalization of the exchanged constant in the Ising 
model, namely Jpotts = '2'Jising-) In the limit of an infinitely long strip {N oo), 
Eq. ( p2[ ) transforms to the duality condition for the largest eigenvalue: 

x[^\k*) = {sinh2K)-^x[^\K). (24) 
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It follows from here that the values of the normalized internal energy in dually 
conjugated points {K and K*) are connected by the relation 



ul{K*)^^ul{K)~2u^{K) (25) 

in which the additive term uo (= coth2iir) does not depend on L. Another 
important feature, this time related to isotropy of a lattice, is that the dually 
conjugated points K and K* confluent into one point at criticality: 

K* = K = Kc. (26) 

Using condition ( p3| ) we find that at the critical point Kc = (l/2)ln(l + y/2) 
the derivative dK* / dK\^ = — 1. Consequently 

ul{K,) = uoiK,) = V2. (27) 

Thus, the critical internal energy per site ul{Kc) of an Ising cylinder with 
isotropic square cells obeys the condition ([l^ ) for all L = 1,2,.... That, in 
turn, leads to the validity of Eqs. (|l|) an d (|17| ). 

In an analogous way, Eqs. (^) and (|17|) can be derived for other isotropic 
spin models partition functions which satisfy a functional equation like 

Zl{K*) = [g{K)]^ZL{K). (28) 

In those cases when the model is self-dual but the critical manifold is a line 
or a surface (as, e.g., for anisotropic lattices), Eqs. (|l|) and ( p7| ) no longer hold. 

This is not difficult to prove if we consider again a two-dimensional Ising 
model. For the anisotropic square lattice the duality condition now looks as 

Xi'-^Kl.Kl) = (sinh2if, sin\i2Ky)-^'^\[''\K^,Ky) (29) 

with 

tanhif* = e"2if„ ^nd tanhi^^ = e^^^:.^ (30) 
From this it follows that on the critical line 



sinh 2K^ sinh 2Ky = 1 (31) 

a condition ( |29| ) relates the values of the free energy in distinct (dually conju- 
gated) points (Kx,Ky) and {Ky,Kx): 

hiK,, Ky) = fL{Ky,K,) + i ln(sinh 2K, sinh 2Ky). (32) 

This circumstance prevents fulfillment of the validity of Eqs. (|l|) and (|^ which 
identify the internal energies at the same point. 
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The critical internal energy density of a strip Lx oo cut out from anisotropic 
lattice depends on the size L. This is easy to verify if, using the results of the 
section 2, one calculates the values ui{Kc) and U2{Kc) for the anisotropic Ising 
and Potts lattices. 

On the other hand, we can establish the same property if we take the Onsager 



solution 1 24 for the two-dimensional Ising model. The dominant eigenvalue of 
a transfer matrix of cylinder L x oo with spatially anisotropic interactions is 
equal to 

x[''\K; a) = (2sinh2X)^/2gxp[(7i + 73 + . . . + 72L-i)/2], (33) 

where a = Jy / Jx 'is the lattice anisotropy parameter and 7^ are positive solutions 
of the equations 

, „ , sinh2aiir j nrx 

cosh jr = cosh 2aK coth 2K . , cos — . (34) 

' sinh 2K \L J ^ ^ 

From this, for the internal energy per site, we obtain 

The dependences ^r{K) have a smooth extremum (minimum) which lies in the 
isotropic case (a = 1) exactly a,t K = and therefore 



dK 



= [r^ 0). (36) 



As a result, the second term in Eq. (^5|) disappears and the critical internal 
energy ceases to depend on L. When a 7^ 1 then Eq. (|3^) does not have a place 
and UL{K]a) depends on the strip width in a complicated way. This explains 
the failures of exact calculations of Kc from Eq. (|l|) in the anisotropic Ising 
lattice jl^. 

Ending the section we note that, in spite of Eqs. (^ and (p7|), Eq. ( p^ 
cannot be deduced from a dual invariance of system. 



5 Conclusions 

Using a group-theoretical approach we obtained the exact analytical solution for 
the double Potts chain with the Hamiltonian (|^) . The solution gives a possibility 
to examine Eq. (^) for a large number of models both with Ising {q — 2) and 
arbitrary Potts spins (including non- integer q) . The validity of Eq. (^ for the 
triangular Potts lattice with pure three-site interactions in alternate triangular 
faces was established. 
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In this paper it was also shown that Eqs. and ( ]17| ) are a consequence of 
a dual symmetry of models for which the critical point coincides with its dual 
image. 

As far as the author knows, the inverse theorem has not been proved. Du- 
ality plus isotropy or, more correctly, self-conjugation of a critical point are not 
necessary conditions for Eq. (|^). Therefore, generally speaking, there can exist 
systems which are not invariant under dual transformation or a combination 
of dual and star-triangular ones, but for which all amplitudes of finite-size cor- 
rections to the critical internal energy (or to some other quantity) are equal to 
zero. 
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Figure captions 



Fig. 1. Geometry of two-site couplings in the double g-state Potts chain with 
Sq symmetry. 

Fig. 2. Fragment of Potts lattice with three-site interactions in alternate 
triangular faces (shaded). 
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